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SUMMARY 


This  report  presents  an  experimental  study  of  the  stability  of  thin-walled 
unstiffened  circular  cylindrical  shells  under  nonuniformly  distributed  axial 
load.  It  demonstrates  that  such  shells  will  buckle  when  the  stress  at  any 
point  on  their  surface  reaches  the  critical  value  for  uniform  axial  compres¬ 
sion. 

It  is  furthermore  demonstrated  that  statistical  data  with  regard  to  Bhell 
buckling  can  be  obtained  from  limited  tests  on  single  specimens. 
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INTRODUCTION 


The  behavior  of  circular  cylindrical  shells  under  uniformly  distributed 
axial  compression  is  a  subject  which  has  received  much  attention  from 
theoreticians  and  experimentalists  alike.  However,  in  practice  many 
cylindrical  shells  are  used  under  conditions  in  which  this  idealized  load 
distribution  is  neither  achieved  nor  even  approached.  Many  cases  can  be 
cited  in  real  engineering  structures  where  concentrated  loads  are  diffused 
into  a  shell,  and  many  other  cases  in  which  other  nonuniform  systems  are 
encountered  can  be  found  readily.  An  important  group,  for  example,  is 
that  in  which  axial  compression  is  associated  with  flexure. 

In  1932,  Flugge"*-  investigated  the  buckling  of  cylindrical  shells  under  com¬ 
bined  bending  and  compression  and  derived  the  interaction  relationship  for 
a  particular  radius-to-thickness  ratio  and  a  special  longitudinal  buckle 
half  wavelength  radius  ratio.  For  the  case  which  he  considered,  he  demon¬ 
strated  that  the  ratio  of  the  maximum  critical  stress  for  bending  alone  would, 
within  the  terms  of  his  analysis,  be  33  percent  greater  than  the  critical 
stress  for  pure  compression.  Timoshenko,^  in  his  theory  of  elastic  stab¬ 
ility,  referenced  this  work  but  omitted  the  required  qualification  with 
regard  to  buckle  wavelength.  In  193^  Donnell^  presented  his  paper,  "A  New 
Theory  for  the  Buckling  of  Thin  Cylinders  Under  Axial  Compression  and  Bending", 
to  the  Fourth  International  Congress  for  Applied  Mechanics.  He  noted,  as  a 
result  of  this  work,  that  in  pure  bending,  buckling  takes  place  when  the 
stress  at  a  point  in  the  cylinder  wall  about  degrees  to  the  neutral  axis 
rises  to  the  value  which  produces  failure  in  uniformly  stressed  specimens. 
Despite  this  conclusion,  the  validity  of  Flugge's  specific  result  has  been 
accepted  by  many  until  quite  recently. 

In  1959,  Abir  and  Nardo^  published  their  paper  on  "Thermal  Buckling  of 
Circular  Cylindrical  Shells  Under  Circumferential  Temperature  Gradients". 

They  reached  the  conclusion  that  the  axial  buckling  stress  under  variable 
thermal  stress  conditions  is  close  to  the  critical  stress  of  the  cylinder 
when  it  is  subjected  to  uniform  axial  compression,  if  the  variation  of  the 
intensity  of  the  thermal  stress  is  not  large  within  a  half  wavelength  of 
the  buckling  pattern.  This  work  was  followed  by  a  study  made  by  Bijlaard 
and  Gallagher*  on  the  "Elastic  Stability  of  a  Cylindrical  Shell  Under 
Arbitrary  Circumferential  Variation  of  Axial  Stress".  The  prime  conclusion 
of  this  research  was  that  the  cylindrical  shell  buckles  when  the  stress  at 
some  point  reaches  the  critical  stress  for  a  cylindrical  shell  under  uniformly 
distributed  axial  compression.  Seide  and  Weingarten,  in  1961,  reported 
their  work,  "The  Buckling  of  Circular  Cylindrical  Shells  Under  Pure  Bending." 
Their  analysis  showed  that  the  maximum  critical  bending  stress  would,  for 
all  practical  purposes,  equal  the  critical  compressive  stress. 

Experimental  work  to  support  these  theoretical  considerations  has  been 
remarkably  lacking.  It  is  true  that  in  1933  Lundquist'  published  a  note  on 
"Stress  Tests  of  Thin-Walled  Duraluminum  Cylinders  in  Pure  Bending",  while 
in  more  recent  times  Suer,  Harris,  Skene,  and  Benjamin0  published  a  partially 
experimental,  partially  theoretical  paper  which  treated  the  bending  stability 
of  thin-walled  unstiffened  circular  cylinders,  including  the  effects  of 
internal  pressure.  The  most  recent  contribution  is  an  experimental  study 
made  by  Heise .9  This  work,  which  was  published  after  the  present  investi- 
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gation  was  well  In  hand,  bears  direct  comparison  with  part  of  the  work 
reported  herein.  Heise  uses  a  very  similar  circular  traverse  system  in 
the  determination  of  the  quality  of  his  test  vehicles.  The  main  difference 
lies  in  the  definition  of  the  buckling  load.  However,  there  still  appears 
to  be  a  great  need  for  experimental  studies  to  provide  practical  confirmation 
of  existing  theoretical  opinions. 

In  view  of  the  practical  application  of  such  information,  a  program  of 
experimental  studies  was  undertaken  in  this  field.  Investigations  were  made 
not  only  of  those  cases  in  which  there  is  a  smooth  variation  in  load  distri¬ 
bution  but  also  of  those  in  which  there  are  sharp  discontinuities  in  the 
loading  actions.  Studies  were  made  of  load  distributions  in  which  flexure 
as  well  as  compression  was  present. 

Recent  researches  on  the  behavior  of  cylindrical  shells  under  uniform  axial 
compression  have  shown  in  a  positive  fashion  that  the  experimental  approaches 
of  the  past  have  been  far  too  restrictive  in  outlook.  New  and  powerful 
techniques  have  been  and  are  being  developed  to  inprove  this  situation.  A 
new  philosophy  with  regard  to  experimental  studies  on  the  buckling  of  shell 
bodies  has  recently  been  presented  by  Horton.  The  work  which  he  reports 
has  shown  that  there  are  two  basic  statistical  approaches  which  can  be  taken 
in  such  studies:  the  first  is  to  consider  a  population  of  nominally  identical 
specimens,  and  the  second  is  to  consider  a  population  of  buckles  within  the 
test  vehicle  itself. 

In  the  work  reported  here,  both  approaches  have  been  made.  They  lead  to 
the  same  conclusion;  namely,  that  irrespective  of  the  nature  of  the  distri¬ 
bution  of  load  around  the  periphery  of  the  shell,  buckling  takes  place  when 
the  stress  reaches  the  level  which  would  be  critical  for  uniform  load  con¬ 
ditions. 

With  regard  to  the  classic  approach  (namely,  that  of  a  population  of  nomin¬ 
ally  identical  specimens),  it  was  reasoned  that  if  the  tests  were  performed 
using  shell  specimens  of  consistently  high  quality,  both  in  uniformity  of 
geometry  and  in  material  properties,  then  the  initial  imperfections  would 
most  likely  be  randomly  distributed  and  the  buckling  load  resulting  from  the 
test  of  a  single  specimen  could  be  regarded  as  one  member  of  a  population 
of  buckling  loads  for  that  test  configuration.  It  was  further  reasoned  that 
statistical  properties  of  the  population,  such  as  mean  value  and  distribution 
of  values  about  the  mean,  could  be  estimated  from  a  random  sample  of  buckling 
loads  obtained  from  identical  tests  of  several  specimens.  The  effect  on 
shell  buckling  strength  of  two  different  distributions  of  loading  would  then 
be  detected  from  a  comparison  of  the  statistical  properties  of  their  respec¬ 
tive  populations  of  buckling  loads. 

It  is  reasonable  to  expect  that  for  shell  specimens  such  as  described  above, 
the  buckling  loads  for  a  given  test  configuration  would  be  distributed 
normally;  i.e.,  according  to  the  normal  lay  of  error.  In  his  treatment  of 
the  subject  of  scientific  research,  Wilson11  points  out  that  it  appears  to 
be  safe  to  use  the  normal  law  for  observations  where  four  or  more  sources 
of  error  enter  with  about  equal  weight.  This  would  be  the  expected  case 
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for  the  reported  buckling  tests.  For  cylinders  manufactured  to  high  and 
carefully  controlled  quality  standards,  errors  introduced  by  initial  imper¬ 
fections  in  the  shell  body  are  expected  to  be  of  equal  weight  with  those 
errors  introduced  through  the  testing  machine,  through  the  loading  heads, 
and  through  repeated  adjustment  of  the  overall  setup  between  the  tests  in  a 
given  sequence. 

The  mathematical  statement  of  the  normal  law  is  well-known  and  for  application 
in  the  present  instance  will  be  written  as 

-(P-P)2 

M  =  _JL_e  2a2  *p  (l) 

N  o/2IT 


where  AN/N  is  the  fraction  of  population  of  observed  buckling  loads  having 
values  in  the  load  interval  P  +  AP,  and  where  P  and  a  are  the  population 
mean  buckling  load  and  standard  deviation,  respectively.  Figure  1  shows 
graphically  the  statistical  properties  of  a  population  of  buckling  loads 
distributed  in  a  normal  manner.  Thus,  for  such  a  population  it  appears 
that  a  comparison  of  the  P  and  a  values  provides  a  sufficient  basis  for 
showing  the  effect  of  load  distribution  on  the  buckling  strength  of  a  cylin¬ 
drical  shell  structure. 

The  problem,  of  course,  from  an  experimental  point  of  view  is  the  manufacture 
and  acquisition  of  a  large  number  of  test  vehicles.  This  is  a  question  of 
such  importance  that  it  merited  careful  investigation.  A  novel  solution  was 
found. in  the  can  manufacturing  industry. 

The  manufacture  of  cans  is  a  fully  automated  process.  It  has  been  developed 
to  the  extent  that  such  devices  are  readily  and  economically  obtainable. 

They  are  vehicles  which,  to  a  large  degree,  can  be  considered  of  a  similar 
character.  There  is,  in  fact,  an  almost  infinite  population  available,  if 
required . 

There  is,  as  we  mentioned  before,  another  approach  to  this  problem:  the 
approach  which  considers  a  population  of  buckles  within  the  test  vehicle 
itself.  It  is  to  this  approach  that  we  now  direct  our  attention.  Buckle 
population  studies  can  be  made  only  if  it  is  possible  to  generate  a  population 
of  sufficient  size  and  character  to  be  studied  statistically.  The  photograph 
given  in  Figure  2  shows  a  fully  developed  buckle  population  in  a  thin-walled 
cylindrical  shell  under  uniform  axial  compression.  As  you  will  see,  this  is 
a  couple te  buckling  of  the  shell.  The  buckles  are  all  identical  in  character. 
In  Figures  3  through  5>  various  stages  of  a  test  are  carried  out  with  this 
process;  the  resulting  load  population  curve  shown  in  Figure  6  is  recognized 
to  be  the  normal  logistic  curve.  It  is  essential  in  the  operation  of  this 
kind  of  test  procedure  that  precautions  are  taken  to  ensure  that  the  buckle 
process  remains  elastic  throughout  the  test.  This  leads  to  the  requirement 
that  the  depth  to  which  a  buckle  is  permitted  to  develop  be  restricted. 
Restriction  is  acconplished  by  means  of  an  inner  mandrel.  A  cross  section  of 
the  test  vehicle,  together  with  this  restraining  mandrel,  is  given  in  Figure  7* 
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The  developments  referred  to  in  the  previous  paragraph  are  developments 
which  apply  only  to  the  case  of  pure  axial  conpression.  However,  in  this 
work  we  are  concerned  with  nonuniform  distributions  of  axial  load,  and  the 
question  which  faces  us  is.  Can  this  technique  be  modified  to  deal  with 
this  problem?  When  we  consider  the  conplexity  of  obtaining  nonuniform 
distributions  by  means  of  a  series  of  individual  loads,  or  by  loading 
only  a  snail  portion  of  the  perimeter,  we  find  that  if  we  wish  to  attempt 
to  run  a  large  number  of  tests  on  a  single  specimen,  the  design  of  the  neces¬ 
sary  test  gear  is  expensive.  However,  we  recognize  that  nonuniform  distri¬ 
butions  of  load  can  be  readily  and  sinply  obtained  by  using  a  combination 
of  compression  and  flexure.  This  makes  an  extremely  simple  rig  and  makes  it 
possible  to  study  the  effect  of  quality  on  buckle  characteristics  together 
with  the  effect  of  nonuniform  distributions  of  load  on  the  overall  buckle 
behavior.  Provided  a  specimen  can  be  designed  in  such  a  manner  that  its 
properties  are  not  influenced  by  repeated  tests,  we  should  be  able  to  acquire 
a  large  amount  of  statistical  evidence  from  a  single  specimen.  The  method 
of  acconplishing  this  and  the  details  of  the  results  obtained  are  given  in 
the  main  body  of  the  report. 
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Figure  1.  A  Normal  Distribution  of  Buckling  Loads 
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Experimental  Values  are  Given  in  Table  VII 


Figure  6.  Observation  of  Buckle  Formation  With  Increasing  Load 
for  the  Cylinder  of  Figure  2  With  Central  Loading. 
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Figure  7«  Cross  Section  of  Test  Vehicle 

Together  With  Restraining  Mandrel. 
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SHELL  SPECIMEN  FOR  FIRST  SERIES  OF  TESTS 


In  the  first  series,  3^9  cylindrical  shell  specimens  were  tested.  These 
were  all  manufactured  by  mass-production  machinery  so  that  variations  in 
geometry  from  specimen  to  specimen  would  be  confined  within  the  close 
toleiance  limits  essential  for  successful  operation  of  such  equipment. 

Blanks  for  the  cylindrical  bodies  were  punched  from  twice  cold-rolled  steel 
sheet  supplied  from  the  mill  in  the  form  of  rolled  strips.  In  one  continuous 
operation  these  blanks  were  formed  into  a  cylindrical  shape,  the  edges  were 
prepared  and  joined  by  soldering,  and  a  sizing  operation  on  the  completed 
cylinder  was  performed.  Edge  restraint  was  provided  through  a  subsequent 
operation  in  which  the  ends  of  the  cylinders  were  joined  with  a  steel  sheet 
metal  cover  by  the  same  process  as  is  used  for  beverage  cans.  The  shell 
specimens  were  then  ready  for  testing. 

Completed  cylinders  were  2.63  inches  in  diameter  and  4.75  inches  in  length. 
Shell  thickness  was  .0058  inch,  giving  an  R/t  ratio  of  226.  A  typical 
stress-strain  curve  for  the  steel  shell  material  is  given  in  Figure  8. 
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Figure  8.  Stress-Strain  Curve  for  Material  of  the  Series  2  Shell  Specimens. 


TEST  PROCEDURE  FOR  FIRST  SERIES  OF  TESTS 


The  general  arrangement  for  testing  is  Bhovm  in  Figure  9*  Loading  was 
applied  by  either  a  standard  Baldwin-Lima-Hamilton  or  a  Tinius-Olsen  testing 
machine,  both  of  which  are  of  60,  OOO-lb  capacity.  The  machine  load  was 
applied  through  a  spherical  loading  block  in  order  to  minimize  the  effect 
of  slight  variations  in  parallelism  of  the  specimen  ends. 

As  mentioned  previously,  three  general,  cases  of  loading  were  used.  These 
are  shown  schematically  in  Figures  10  through  14.  For  each  case  a  special  set 
of  loading  plates  was  used  to  distribute  the  test  load  around  the  shell 
perimeter.  These  plates  were  made  of  l/4- inch-thick  structural  grade  aluminum 
and  during  a  test  they  contacted  the  specimen  along  the  rings  formed  by 
the  fabrication  Joint  between  the  cylinder  and  the  end  covers.  For  exanple, 
the  set  of  loading  plates  for  the  case  of  three-place  symmetrical  loading 
(Figures  12  and  14)  consisted  of  five  pairs  of  plates,  each  with  three  lobes 
centered  120  degrees  apart.  One  pair  each  was  provided  for  distributing 
the  load  along  18-,  40-,  60-,  80-,  and  100-degree  segments  of  the  end  perim¬ 
eters.  Loading  was  applied  symmetrically  with  respect  to  the  midplane 
normal  to  the  shell  axis.  The  manner  of  arranging  a  set  of  plates  to 
give  a  desired  distribution  of  load  is  clear  from  Figure  9* 

The  end  rims  mentioned  above  possessed  a  degree  of  flexural  stiffness  and 
could  resist  bending  out  of  the  end  plane .  Advantage  was  taken  of  this  fact 
to  create  two  variations  on  each  of  the  cases  of  symmetrically  distributed 
loading.  These  are  depicted  in  Figures  11  and  13  and  Figures  12  and  14  f or 
the  cases  of  two-place  and  three-place  symmetrical  loading,  respectively. 

The  distribution  labeled  Type  B  occurs  when  the  end  rims  are  continuous; 
that  is ,  when  the  flexural  stiffness  of  the  rim  participates  in  the  transfer 
of  load  from  loading  plate  to  shell.  By  the  simple  expedient  of  cutting 
through  the  rim  at  the  ends  of  a  loaded  segment,  the  distribution  of  Figure 
10,  labeled  Type  A,  was  caused  to  occur.  This  cutting  was  done  with  a  thin, 
high-speed,  abrasive  wheel  without  affecting  the  shell  body. 

That  the  Type  A  and  Type  B  variations  of  load  distribution  did  occur  was 
shown  with  the  aid  of  SR- 4  electrical  resistance  strain  gages.  A  number  of 
these  having  l/8-inch  grid  length  were  placed  around  an  end  of  a  shell 
specimen  so  that  peripheral  distribution  of  longitudinal  strain  could  be 
measured.  Different  load  plates  were  used  to  distribute  loading  through  both 
a  continuous  and  a  cut  rim.  Strain  measurements  from  these  tests  provided 
the  data  from  which  the  load  distributions  shown  in  Figures  10  through  14 
were  constructed. 

After  a  specimen  had  been  set  up  for  tests,  the  machine  load  was  gradually 
increased.  Typically,  the  load  would  be  observed  to  rise  smoothly  until 
failure  by  buckling  occurred,  after  which  it  would  suddenly  fall  off  and 
stabilize  at  a  lower  value.  If  the  machine  loading  process  was  continued, 
the  load  would  rise  slowly  until  it  was  slightly  above  this  initial  stable 
falloff  value,  and  then  a  second  buckling  and  load  falloff  sequence  would 
occur.  As  testing  progressed,  the  initial  falloff  came  to  be  considered 
a  significant  characteristic  of  the  buckling  failure,  data  on  which  might 
be  of  velue  in  planned  future  experiments. 
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The  data  recorded  for  each  test  included  the  buckling  load,  the  location 
of  the  failure,  the  sound  of  the  buckle  formation  (whether  a  sharp  or  a 
soft  snap),  and,  in  most  instances,  the  falloff  load.  Tables  I  through  V 
contain  all  recorded  test  data  from  the  shell  specimens  used  in  the 
experimental  program.  Reference  is  made  in  each  table  to  the  corresponding 
load  distribution  diagram  of  Figures  11  through  14.  Typical  buckling 
failures  are  shown  in  Figure  15 . 
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Figure  15.  Typical  Buckling  Failures,  Series  1  Shell 
Specimens. 


RESULTS  OF  THE  FIRST  SERIET 


The  critical  forces  for  unsymmetrical  loading  of  the  type  designated  A  in 
the  test  are  presented  in  Table  I.  The  values  for  tvo-place  symmetrical 
loading  under  the  same  condition  of  edge  fixity  are  listed  in  Table  II,  while 
the  corresponding  levels  for  three-place  loading  are  given  in  Table  III. 
Probability  plots  and  cumulative  distribution  curves  for  two  representative 
samples  of  the  data  outlined  above  sure  shown  in  Figures  16  and  17» 

The  data  relative  to  tvo-place  symmetrical  and  three-place  symmetrical 
loading  with  edge  fixity  Type  B  are  shown  in  Tables  IV  and  V,  respectively. 

A  probability  plot  and  cumulative  distribution  curve  for  an  appropriate 
sample  of  this  information  is  exhibited  in  Figure  18.  The  mean  values  of 
compressive  load  to  produce  instability  are  plotted  against  the  percentage 
of  perimeter  loaded  for  types  A  and  B  edge  conditions  in  Figures  19  and  20. 
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Type  A,  Loaded  Fraction  of  Perimeter  *  .50. 


Figure  17.  Results  of  Testing  a  Sample  of  25  Cylindrical  Shells,  Unsymaetrical  Distributed 
Loading  -  Type  A,  Loaded  Fraction  of  Perimeter  *  .86. 
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Figure  18.  Results  of  Testing  a  Sample  of  50  Cylindrical  Shells,  Two-Place  Symmetrical 
Loading  -  Type  B,  Loaded  Fraction  of  Perimeter  =  .50. 


SAMPLE  MEAN  BUCKLING  LOAD,  LB 


BASIC  DATA  CONTAINED  IN  TABLES,I«1I  ,•  ID 


0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

FRACTION  OF  PERIMETER  LOADED 


Figure  19.  Buckling  Load  as  a  Function  of  Fraction  of  Perimeter 
Loaded  for  Series  1  Shell  Specimen  -  Type  A  Loading. 
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SAMPLE  MEAN  BUCKLING  LOAD,  LB 


BASIC  DATA  CONTAINED  IN  TABLES  E7  AND  7 


Figure  20,  Buckling  Load  as  a  Function  of  Fraction  of  Perimeter 
Loaded  for  Series  1  Shell  Specimen  -  Type  B  Loading. 
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TABLE  I.  UNSY METRICAL  DISTRIBUTED  LOADING  -  TYPE  A 
Testing  Machine:  60, OOO-lb  Tinius-Olsen 
(See  Figure  10) 


Chord  Perimeter  Buckle  Buckle  Data  Falloff 


Shell  No. 

(in) 

(pet)  Load  (lb) 

1  Location 

Snap 

Soft 

Load 

GROUP  1  (10  Tests) 

S4lU 

2-1/4 

.33  630 

Top 

X 

440 

S42U 

510 

Bottom 

X 

340 

S43U 

590 

Top/Middle 

X 

- 

S44u* 

620 

Upper  l/3 

X 

460 

S45U 

740 

Top 

X 

340 

S46U 

580 

Top 

X 

300 

S47U 

750 

Top 

X 

380 

S48u 

600 

Top/Middle 

XX 

480/320 

S49U 

540 

Top 

X 

420 

S50U 

640 

Top 

X 

320 

GROUP  2  (40) 

S1U 

2-5/8 

.50  840 

Top 

XX 

720/580 

S2U 

990 

Top /Twiddle 

XX 

820/640 

S3U* 

900 

Bottom/Middle 

XX 

840/680 

s4u 

930 

Top 

X 

640 

S5U 

900 

Top 

X 

700 

S6u 

850 

Middle 

XX 

770/520 

S7U 

840 

Top/Middle 

XX 

780/520 

S8u 

810 

Top/Middle 

X 

560 

S9U 

820 

Top 

X 

- 

S10U* 

740 

Top 

X 

500 

snu* 

1080 

Middle 

XX 

600/460 

S12U 

1000 

Top 

XX 

820/500 

si3u 

760 

Top 

X 

460 

Sl4u 

890 

Top 

X 

520 

S15U 

1000 

Top 

X 

520 

si6u 

810 

Top 

X 

480 

S17U 

690 

Top 

X 

620 

Sl8u 

840 

Bottom 

XX 

680/480 

S19U 

700 

Top 

X 

470 

S20U 

790 

Top 

X 

540 

S21U 

880 

Top 

X 

430 

S22U 

710 

Middle 

X 

600 

S23U 

930 

Top 

X 

46c 

S24U 

810 

Top 

X 

540 

S25U 

930 

Top 

X 

460 

S26U 

1010 

Bottom 

X 

540 

S27U 

940 

Top 

X 

420 

S28U 

800 

Top/Middle 

XXX 

750/750/680 

S29U* 

920 

Bottom/Middle 

X 

520 

S30U 

700 

Top /Middle 

X 

540 
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Shell  No. 


TABLb  I  -  continued 


Chord  Perimeter  Buckle  _ Buckle  Data _  Falloff 

(in)  (pet)  Load(lb)  Location  Snap  Soft  Load 


GROUP  2  (Cont'd) 


S31U 

S32U 

S33U 

S34U 

S35U 

S36U 

S37U* 

S38U 

S39U 

s4ou 


860  Top/Middle 

810  Top 

760  Top/Middle 

910  Top 

890  Top 

900  Middle 


GROUP 


2-1/4 


500 

520 

620 

680 

600 

820/720/5 


99  0 

Bottom 

X 

500 

94o 

Top 

X 

460 

990 

Top 

X 

520 

990 

Top 

X 

440 

’  3  do! 

1120 

L 

Top/Middle 

X 

X 

780 

1020 

Top 

X 

680 

1310 

Top 

X 

580 

1100 

Top 

X 

540 

1390 

Top 

X 

600 

1210 

Top 

X 

630 

mo 

Top 

X 

740 

1260 

Top 

X 

X 

780 

1220 

Top/Middle 

X 

mo 

1050 

Top/Middle 

X 

- 

GROUP  4  (25 


S6lU 

S62U 

S63U 

s64u* 

S65U 

s66u 

s6tu 

s68u* 

S69U 

S70U 

S71U 

S72U 

S73U* 

S74u 

S75U 

S76U 

S77U 

S78U 

S79U 

S80U 


1-1/8 


1410 

Top /Middle 

X 

900 

1680 

Top/Middle 

X 

830 

1540 

Top  /toddle 

X 

900 

1800 

Top /Middle 

X 

1020 

1790 

Top 

X 

760 

1730 

Top 

X 

660 

1560 

Top 

X 

64o 

1690 

Top/Middle 

X 

760 

1710 

Top 

X 

740 

1620 

Middle 

X 

1060 

1360 

Top 

X 

800 

1500 

Top /Middle 

X 

960 

1590 

Middle 

X 

960 

1470 

Top 

X 

780 

1460 

Middle 

X 

1280 

1300 

Top /Middle 

X 

800 

1630 

Top 

X 

820 

1450 

Top 

X 

740 

1200 

Top 

X 

740 

ltoo 

Top/Middle 

X 

840 
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TABU;  I  -  continued 

„  ,,  ,,  Chord 

Shell  No •  ^ 

Perimeter  Buckle  Buckle  Data 

(pet)  Load  (lb)  Location  Snap  Soft 

Falloff 

Load 

GROUP  k  (Contact) 


S81U 

1510 

Top/Middle 

X 

780 

S82U 

1500 

Top 

X 

7to 

S83U 

1520 

Top 

X 

760 

s8to 

15  to 

Top 

X 

700 

S85U 

1270 

Top /Middle 

X 

7to 

GROUP  5  (2) 

S86u 

0 

0 

160 

X 

S37U 

0 

0 

180 

X 

- 

*  Excellent  Buckle  Pattern 
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TABLE  II.  TWO-PLACE  SYMMETRICAL  DISTRIBUTED  LOADING  -  TYPE  A 
Test  Machine:  60, 000-lb  Tlnius-Olsen 
(See  Figure  11 ) _ _ _ _ 

Buckle  Buckle  Data 

Shell  No.  Width  Load  Location  Snap  Soft  Falloff  Load 

GROUP  1  (10  Tests) 


S103a 

1 

580 

Top 

X 

200 

S104a 

580 

Top /Bottom 

X 

- 

S105a 

530 

Top 

X 

320 

S106a 

470 

Top 

X 

340 

S107a 

500 

Top 

X 

380 

S108a 

490 

Top 

X 

390 

S109a 

480 

Top 

X 

360 

S70 

600 

Top 

X 

- 

S71 

480 

Top /Middle 

X 

- 

S72 

500 

Top/Middle 

X 

• 

GROUP  2  (50) 

S80a 

1-7/8 

1020 

Top 

X 

580 

S8la 

1100 

Middle 

X 

720 

S82a 

960 

Top 

X 

500 

S83a 

1040 

Bottom 

X 

730/520 

S84a 

980 

Top 

X 

520 

S85a 

1040 

Top 

X 

520 

S86a 

1020 

Bottom 

X 

790 

S87a 

1020 

Top 

X 

560 

S88a 

810 

Top/Middle 

X 

X 

660 

S89a 

990 

Top/Middle 

X 

700 

S90a 

820 

Top 

X 

660 

S91a 

1020 

Top 

X 

540 

S92a 

1010 

Top 

X 

600 

S93a 

840 

Top 

X 

500 

S94a 

1010 

Top 

X 

570 

S95a 

910 

Middle 

X 

800 

S96a 

920 

Bottom/Middle 

X 

780 

S97a 

800 

Middle 

X 

760 

S98a 

1100 

Top 

X 

660 

S99a 

980 

Top 

X 

X 

600 

SlOOa 

950 

Middle 

X 

860/710 

SlOla 

1010 

Top /Middle 

XX 

960/670 

S102a* 

990 

Middle 

690 

S80b 

1080 

Top 

X 

550 

S8lb 

1120 

Top/Middle 

X 

640 

S82b 

1070 

Top 

X 

580 

S83b 

950 

Top 

X 

600 

S84b 

1090 

Top 

X 

600 
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TABI.b  II  -  continued 

—  '  Buckle  Buckle  Data 

Shell  Wo.  Width  Load  Location  Snap 

GROUP  2  (Cont'd) 


S85b 

1220 

Top 

X 

S86b 

980 

Middle 

X 

S87b 

1110 

Top 

X 

S88b 

1130 

Top 

X 

S89b 

1030 

Top 

S90b 

1180 

Top 

X 

S91b 

1170 

Top/Middle 

X 

S92b* 

1220 

Top/Middle 

X 

S93t 

860 

Top/Middle 

X 

S94b 

1120 

Bottom 

X 

S95b 

970 

Top 

X 

S96b 

890 

Top 

S97b 

1120 

Top/Middle 

X 

S98b 

1080 

Top 

X 

S99b 

910 

Bottom/Middle 

X 

SlOOb 

1090 

Top 

X 

SlOlb 

1150 

Top 

X 

S102b 

1090 

Top 

X 

S103b* 

1050 

Bottora/Middle 

X 

S104b 

1170 

Top 

X 

S73 

980 

Top 

X 

S74 

1010 

Top 

GROUP  3  (10) 

S70a 

2-7/16  1310 

Top 

X 

S71a 

1570 

Bottom 

X 

S72a 

1420 

Top 

X 

S73a 

1440 

Top 

X 

S74a 

1390 

Top 

X 

S75a 

1400 

Bottom 

X 

S76a 

1380 

Top /Middle 

X 

S77a 

1420 

Top 

X 

S78a 

1470 

Top 

X 

S79a 

1420 

Top 

X 

*  Excellent  Buckle  Pattern 
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Soft  Falloff  Load 


580 

800 

640 

x  660 

600 
640 
700 
680 
620 
900 
840 

x  600 

700 
620 
740 
580 
620 
600 
900 
600 


x 


1120 

780 

760 

780 

1120 

980/850 

1200/840/720 

780 

900 
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TABLIi  III  -  continued 


Shell  No. 

Buckle 
Cl  Load 

Buckle  Data 
Location  Snap 

USB 

Falloff  Load 

GROUP  3  (Cont'd) 

S62c 

1610 

Top /Middle 

X 

nto 

S63c 

1580 

Top 

X 

1100 

S64c 

1700 

Top 

X 

920 

S65c 

15  ^ 

Bottom/Middle  x 

1360 

S66c 

1690 

Top 

X 

1120 

s67c 

1800 

Top /Middle 

X 

1060 

S68c 

1260 

Bottora/Middl 

e  x 

X 

- 

S69c 

1630 

Middle 

X 

1380 

S70c 

1500 

Top/Middle 

X 

1000 

S71c 

1680 

Bottora/Middl 

e  xx 

1500/1260 

S72c 

1600 

Top 

X 

- 

S73c 

1570 

Top 

X 

1150 

S7to 

1600 

Top 

X 

1180 

S75c 

1760 

Top 

X 

860 

S76c 

1660 

Top 

X 

1200 

S77c 

i6to 

Top/Middle 

X 

1210 

S78c 

1680 

Top 

X 

1200 

S79c 

1620 

Top 

X 

1020 

S80c 

1890 

Top/Middle 

X 

1160 

GROUP  4  (3) 

S2U 

80°  1290 

Top 

X 

- 

S25 

1120 

Top 

X 

- 

S26 

1220 

Top 

X 

- 

*  Excellent 

Buckle  Pattern 
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TABLE  IV. 

TWO-PLACE  SYMMETRICAL  DISTRIBUTED  LOADING  - 
Testing  Machine:  60, OOO-lb  Tinius-Oleen 
(See  Figure  13) 

TYPE  B 

i  i 

Buckle  Buckle  Data 

Shell  No. 

Width  Load  Location  Snap  Soft 

Falloff  Load 

GROUP  1  (3  Tests) 

S56 

1"  685 

Bottom  x 

- 

S57 

745 

Bottom  x 

- 

S58 

800 

Bottom  x 

- 

GROUP  2  (£3) 

100 

1-7/8"  mo 

Top /Bottom 

800 

101 

1000 

Bottom/Middle 

750 

102 

1020 

Top /Middle 

730 

103 

1155 

Top 

800 

104 

1025 

Upper  l/3 

725 

105 

1240 

Top 

930 

106 

940 

Top/Middle 

550 

107 

1125 

Top /Bottom 

690 

108 

1170 

Middle 

720 

109 

1055 

Top 

830 

110 

1080 

Top /Middle 

740 

111 

1110 

Top 

900/810 

112 

950 

Bottom/Middle 

700 

113 

970 

Top 

700 

n4 

1200 

Top /Bottom 

860 

115 

1100 

Top /Low  l/3 

670 

116 

1035 

Top /Bottom 

600 

117 

1240 

Bottom/Middle 

1070 

118 

1105 

Top 

730 

119 

955 

Top /Upper  1/3 

700 

120 

980 

Top 

690 

121 

1155 

Top 

630 

122 

1160 

Top 

800 

123 

1130 

Top 

700 

124 

1105 

Top 

730 

125 

1000 

Top 

700 

126 

1090 

Top 

720 

127 

1065 

Top/Middle 

940 

128 

1045 

Top/Upper  l/3  and  Bottom 

800 

129 

1250 

Top 

740 

130 

970 

Top/Middle  l/3 

550 

131 

1145 

Top 

550 

132 

n4o 

Middle 

740 

133 

970 

Top 

690 

134 

1155 

Top 

650 

135 

1150 

Bottora/Lower  l/3 

550 
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|  TABLE  IV  -  continued 

Buckle 

Buckle  Load 

Shell  No. 

Width 

Load 

Location 

Snap 

Soft 

Falloff  Load 

GROUP  2  (Cont'd) 

136 

1050 

Bottom 

550 

137 

990 

Top /Upper  l/3 

620 

138 

1265 

Top 

820 

139 

915 

Top/Middle  1/3 

650 

l4o 

1070 

Top/Bottom 

790 

l4l 

935 

Top 

680 

142 

1175 

Top 

740 

143 

895 

Top 

700 

144 

980 

Top /Upper  l/3 

760 

145 

1015 

Top 

760 

146 

1085 

Top 

650 

147 

1115 

Top 

720 

146 

1160 

Top 

720 

149 

1050 

Top 

730 

S53 

1110 

Top 

- 

S54 

1165 

Top 

- 

S55 

1140 

Bottom 

- 

GROUP  3  (25) 

150 

2-7/16 

1650 

Top /Bottom 

X 

1200/980 

151 

1710 

Top 

X 

960 

152 

1730 

Top 

xxxx 

1060/980/900/760 

153 

1790 

Top /Middle 

xxxxx 

1100/1080/1010/900/800 

154 

1550 

Top/Middle 

XX 

880/700 

155 

1570 

Top 

X 

1000 

156 

1720 

Middle 

XX 

1260/1000 

157 

1510 

Top 

X 

1100 

158 

1620 

Top 

X 

880 

159 

1560 

Top 

X 

900 

160 

1560 

Top /Middle 

X 

1100 

l6l 

1540 

Top /Middle 

X 

950 

162 

1830 

Top 

X 

1040 

163 

1550 

Top 

X 

900 

164 

1580 

Top 

X 

1000 

165* 

1670 

Top 

X 

900 

166 

1610 

Top 

X 

920 

167 

1620 

Top/Middle 

X 

1200 

168 

1430 

Top 

X 

- 

169 

1510 

Top 

X 

1000 

170 

1430 

Top/Middle 

X 

1020 

171 

1540 

Top 

X 

1000 

S50 

1535 

Top 

X 

- 

S51 

1500 

Top 

X 

- 

S52 

1790 

Top 

X 

- 

*  Excellent  Buckle  Pattern 
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TABLE  V.  THREE-PLACE  SYMMETRICAL  DISTRIBUTED  LOADING  -  TYPE  B 
Testing  Machine:  60, OOO-lb  Baldwin-Lima-Hamilton 
(See  Figure  14) 


Shell  No. 

a 

Buckle 

Load 

Buckle  : 
Location 

Data 

Snap  Soft 

GROUP  1  (4  Tests) 

Sla 

18° 

620 

Top 

X 

S2a 

520 

Top 

X 

S3a 

590 

Middle 

X 

S4a 

610 

Top 

X 

GROUP  2  (11) 

Sl6a 

40° 

970 

Top 

X 

S17a 

800 

Middle 

X 

Sl8a 

1000 

Top  1/3 

X 

S19a 

800 

Bottom/Middle 

X 

S20a 

960 

Top 

X 

S21a 

990 

Top 

X 

S22a 

760 

Bottom/Middle 

X 

S23a 

950 

Top 

X 

SI 

940 

Top 

X 

S2 

1000 

Top 

X 

S3 

1030 

Top 

X 

GROUP  3  (9) 

S5a 

60° 

1020 

Bottom/Middle 

X 

S6a 

1130 

Bottom/Middle 

X 

S7a 

1060 

Top 

X 

S8a 

1080 

Top 

X 

S40a 

I28O 

Top 

X 

S4la 

1130 

Bottom/Middle 

X 

Sk2a 

920 

Middle 

X 

Sk3a 

1230 

Top/Middle 

X 

S44a 

1050 

Middle 

X 

GROUP  4  (11) 


S24a 

80°  150C 

Top  x 

S25a 

1360 

Bottom/Middle  x 

S26a 

1320 

Bottom/Middle  x 

S27a 

1400 

Top  x 

S28a 

1300 

Bottom  x 

S29a 

1260 

Top  x 

TABLE  V  -  continued 


Buckle  Buckle  Data 


Shell  No. 

a 

Load 

Location 

Snap  Soft 

GROUP  5  (7) 

S9a 

100° 

1640 

Top/Middle 

X 

SlOa 

lklO 

Top/Middle 

X 

Slla 

1820 

Middle 

X 

S12a 

1840 

Top /Middle 

X 

S13a 

19l*0 

Top/Middle 

X 

Sl4a 

lW*0 

Top/Middle 

X 

S15a 

1800 

Top/Middle 

X 

GROUP  6  (10) 

S32a 

120° 

1920 

Top 

X 

S33a 

1620 

Top 

X 

S3ka 

1690 

Top 

X 

S35a 

2060 

Top/Middle 

X 

S36a 

2080 

Top 

X 

S37a 

1730 

Top 

X 

S38a 

2070 

Top 

X 

S7 

I960 

Top 

X 

S8 

1660 

Top 

X 

S9 

1760 

Middle 

X 
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DISCUSSION  OF  THE  RESULTS  OF  THE  FIRST  SERIES 


Reference  is  made  to  the  information  contained  in  Figures  16  through  18.  The 
quality  of  the  data  is  seen  to  be  high.  The  standard  deviations  vary  from 
6  percent  at  the  high  load  levels,  which  correspond  to  the  greater  loaded 
areas,  to  12  percent  for  loading  over  a  smaller  fraction  of  the  circumference. 
This  variation  is  well  inside  that  frequently  experienced.  Mossakovskii  and 
Smelyi, 12  for  instance,  quote  15  percent  as  being  not  unreasonable. 

The  data  of  Figure  19,  which  encompass  all  load  distribution  for  the  type  A 
conditions,  are  strong  evidence  that  the  critical  stress  which  causes  insta¬ 
bility  is  Independent  of  the  nature  of  the  end  load  distribution.  Statistical 
treatment  and  normal  correlation  procedures  outlined  in  the  appendix  indicate 
that  the  straight  line  oA  represents  this  material  to  a  confidence  level  of 
95  percent.  The  results  depicted  in  Figure  20  are  totally  consistent  with 
those  referred  to  above. 
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THE  SECOHD  SERIES  OF  TESTS 


In  the  second  series  of  tests,  the  specifications  were  much  broader  than 
those  for  the  first  series.  Previously,  only  the  application  of  classic 
procedures  to  determine  the  influence  of  nonunifonnity  of  stress  conditions 
was  of  any  concern.  Preparations  had  been  made  to  use,  if  necessary,  a  very 
large  number  of  specimens  and  to  accept  the  time-consuming  process  of  setting 
up  and  testing  these  specimens,  whereas  in  this  series  the  use  of  a  single 
specimen  to  develop  a  procedure  for  studying  the  quality  of  the  test  vehicle 
and  at  the  same  time  to  study  the  influence  of  combined  compression  and 
flexure  on  the  buckling  characteristics  of  thin-walled  structures  was  desired. 
Hie  study  was  begun  with  the  observation  taken  from  the  first  series;  namely, 
that  when  a  cylindrical  shell  is  loaded  with  discrete  loads  around  its  cir¬ 
cumference,  the  buckling  which  takes  place  occurs  in  regions  which  are  close¬ 
ly  associated  with  the  loaded  region.  Therefore,  there  mast  be  a  greater 
tendency  to  buckle  in  some  regions  of  the  shell  than  in  others  If  a  cylin¬ 
drical  shell  is  loaded  by  an  offset  load.  Thus,  the  shell  can  be  considered 
In  a  slightly  different  light.  The  regions  which  tend  to  buckle  can  be  con¬ 
sidered  to  be  plates  with  undetermined,  but  nevertheless  repeatable,  loading 
and  support  conditions  for  all  similar  overall  load  application  conditions. 
Hence,  a  single  shell  can,  under  these  conditions,  be  looked  upon  as  a  wide 
population  of  test  vehicles,  each  of  which  bears  somp  connection  with  every 
other  similar  specimen.  Clearly,  the  methods  of  manufacture  of  the  shell  mist 
mean  that  individual  elements,  cut  at  random  around  the  periphery,  would 
differ  from  each  other  with  respect  to  the  quality.  Thus,  it  seems  reason¬ 
able  to  consider  that  a  circular  traverse  would  provide  information  relative 
to  the  quality  of  the  shell  in  regions  which  are  defined  as  lengthwise  strips. 
Experience  with  pure  compressions!  buckling  would  indicate  that  if  the  buckle 
load  is  defined  as  that  load  which  produces  maximum  buckle  generation,  then 
each  strip  would  have  the  same  load;  but  the  kurtosis  of  the  buckle  distri¬ 
bution  curve  would  alter  with  quality,  becoming  more  platykurtic  with  improved 
quality.  In  other  words,  a  good  section  would  tend  to  buckle  uniformly  all 
over  at  one  and  the  same  time. 

Thus,  In  the  first  tests  of  the  second  series,  the  circular  traverse  method 
was  investigated  as  a  technique  for  quality  control.  The  second  tests  of 
the  series  were  made  to  determine  the  influence  of  combined  compression  and 
flexure  on  the  buckling  characteristics  of  thin-walled  circular  cylindrical 
shells.  The  case  when  the  offset  of  load  was  zero,  as  closely  as  could  be 
determined,  was  used  to  check  the  quality  control  procedure. 

TEST  SPECIMENS 


The  specimen  used  in  these  tests  was  manufactured  in  the  laboratory  by  the 
following  procedure.  A  thick -walled  tube  of  appropriate  dimension  and 
specification  was  accurately  bored  and  honed.  It  was  then  carefully  turned 
until  the  wall  thickness  was  of  the  order  of  l/32  inch.  When  this  was  com¬ 
plete,  the  shell  was  shrunk  onto  an  accurate  mandrel  whose  coefficient  of 
linear  expansion  differed  from  that  of  the  shell  material.  The  mandrel- 
shell  combination  was  then  turned  between  centers,  using  a  carbide  tool, 
until  the  final  chosen  thickness  was  approximated.  The  final  dimension  was 
achieved  by  lapping  and  polishing. 


The  shell  so  manufactured  vas  assembled  in  loading  heads  and  was  made 
"encastre"  In  the  loading  heads  with  a  low-melting-point  alloy.  A  buckle 
depth-restricting  mandrel  was  arranged  inside  the  shell,  concentric  with  the 
shell.  A  gap  of  approximately  l/8  inch  was  allowed  between  the  head  of  the 
mandrel  and  the  lower  face  of  the  upper  loading  plate.  A  small  hole  per¬ 
mitted  entrapped  air  to  escape  when  the  shell  volume  was  reduced  by  the 
buckle  formation. 

TEST  PROCEDURE 


The  cylindrical  shell  manufactured  in  accordance  with  the  test  specimens  in 
the  preceding  section  was  tested  in  a  60,000-pound-capacity  Tinius-Olsen 
Universal  test  machine  under  combined  flexure  and  compression.  In  order  to 
make  the  setup  process  as  convenient  as  possible,  a  special  jig  was  devised. 
Details  of  this  are  shown  in  Figure  21.  To  assemble  the  test  specimen  in  the 
machine,  the  following  procedure  was  adopted.  The  special  Jig  was  located 
on  the  loading  platform,  and  the  lower  loading  ball  was  positioned.  Next, 
the  test  specimen  was  aligned  on  the  jig,  being  supported  in  the  unloaded 
state  by  the  three  low  rate  springs  seen  in  Figure  21.  The  top  loading  ball 
was  now  sited  in  the  correct  place,  and  the  machine  head  lowered  until 
contact  was  made  with  the  ball.  By  paying  particular  attention  to  detail, 
and  by  checking  and  rechecking,  it  was  found  possible  to  align  the  upper  and 
lower  loading  balls  vertically  with  great  accuracy. 

For  the  first  tests  made,  the  line  of  action  of  the  applied  load  was  kept  at 
a  constant  distance  A  from  the  axis  of  the  shell,  such  that  A/r  =  3/8. 

Eight  positions,  45  degrees  apart,  were  used.  A  buckle  number  load  history 
was  determined  for  each  loading  station. 

The  next  tests  in  this  series  were  repeats  of  this  first  family  except  that 
A/r  was  changed  to  l/2. 

The  final  tests  were  conducted  using  a  radial  traverse  instead  of  the  cir¬ 
cular  traverse  used  for  the  first  and  second  families  of  the  series.  In  this 
radial  traverse,  seven  values  of  A/r  were  used. 

The  cylinder,  els  set  in  the  machine  for  testing,  is  seen  in  Figure  22. 
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TEST  VEHICLE 
DETAILS  AS  SHOWN 
IN  FIGURE  7 


Figure  21.  Cross  Section  of  Series  2  Test  Vehicle  Kith  Restraining 
Mandrel  and  Testing  Jig. 


DISCUSSION  OF  THE  RESULTS  OF  THE  SECOND  SERIES 


The  mechanical  properties  of  the  material  of  the  test  cylinder  are  given  in 
Figure  23 •  The  results  of  the  first  circular  traverse  with  b/r  =  3/8  are 
given  in  Table  VI  and  are  presented  in  graphical  form  in  Figure  24.  The 
buckle  number  load  histories  displayed  therein  follow  normal  cumulative 
distributions.  The  buckle  load  for  each  section,  defined  as  the  load  for 
maximum  buckle  germination,  is  substantially  constant.  The  variation  which 
does  occur  is  so  slight  that  it  muBt  be  considered  to  be  accidental.  It  is 
seen  that  the  main  difference  among  the  several  data  presented  is  that  at  one 
station  the  standard  deviation  is  very  large  relative  to  the  value  of  this 
parameter  at  all  other  stations.  This  Implies,  then,  that  perhaps  this 
section  should  be  the  weakest  section  of  the  shell  or  at  least  the  section  with 
the  greatest  irregularity. 

The  results  of  the  second  circular  traverse  with  b/r  =  l/2  are  also  presented 
in  Table  VI  and  are  displayed  graphically  in  Figure  25.  The  results  sure  in 
1:1  correspondence  with  those  previously  obtained. 

Figure  26  shows  the  variation  in  buckling  load  for  the  circular  traverse  at 
3/8r  and  l/2r  eccentricities.  The  tendency  toward  a  constant  value  for  any 
radial  position  of  the  load  is  clearly  seen. 

In  Figure  27  the  variation  in  standard  deviation  for  these  two  families  of 
tests  is  given  as  a  function  of  the  angular  position  of  the  test  load.  It 
is  immediately  apparent  that  the  maximum  values  of  a  are  located  in  the 
same  regions. 

The  shell  was  now  tested  with  b/r  =  0,  i.e.,  under  pure  axial  compression. 

The  area  of  first  buckling  under  these  conditions  was  found  to  correspond  to 
part  of  the  region  for  which  the  maximum  value  of  standard  deviation  was 
measured.  The  results  of  testing  under  pure  axial  compression  are  given  in 
Table  VII.  These  are  the  data  used  to  plot  the  example  curve  shown  in 
Figure  6. 

If  the  hypothesis  that  initial  imperfection  determines  the  point  of  initial 
buckle  germination  is  accepted,  then  the  results  given  above  indicate  that: 

1.  Under  identical  load  conditions,  nominally  identical  specimens  have  the 
same  buckle  load,  when  this  load  is  defined  as  the  load  which  causes 
maximum  buckle  generation. 

2.  The  kurtosis  of  the  curve  of  rate  of  change  of  number  of  buckles  as  a 
function  of  load  is  a  measure  of  the  degree  of  imperfection. 

For  the  final  family  of  tests,  the  shell  was  set  such  that  the  defective 
region  as  determined  by  the  above -described  procedure  lay  in  opposite 
quadrant  to  the  radial  traverse  line.  With  the  shell  so  located,  a  radial 
traverse  was  maHe  using  values  of  b/r  from  0  to  3/4,  in  l/8  steps.  The 
buckle  number  load  histories  for  these  offsets  are  given  in  Tables  VII  and 
VIII  and  are  illustrated  graphically  in  Figures  28  and  29.  The  corres¬ 
ponding  buckle  patterns  are  shown  in  Figure  30.  From  these  results  it  is 
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seen  that  in  each  case  the  buckle  number  load  histories  are  normal  cumulative 
distributions.  It  is  observed  that  there  is  no  change  in  the  shape,  size, 
or  location  of  the  buckles  formed  in  the  various  cases,  but  that  the  size 
of  the  population  varies  with  load  location.  The  variation  in  population 
size  as  a  function  of  offset  magnitude  is  portrayed  in  Figure  31*  It  is 
clear  from  this  curve  that  when  A/r  is  greater  than  l/5,  a  very  large  per¬ 
centage  of  the  surface  of  the  shell  is  devoid  of  buckles. 

This  observation  suggests  immediately  the  main  advantage  that  might  derive 
from  the  use  of  a  radial  traverse.  It  seems  reasonable  to  conclude,  bearing 
Saint  Venant '  s  principle  in  mind,  that  the  same  information  might  be  obtained 
from  an  incomplete  shell  as  can  be  obtained  from  a  complete  shell  by  use  of 
the  radial  traverse  procedure.  It  is  pertinent  to  remark  that  if  a  seam  is 
to  be  used,  the  location  for  it  is  clear  and  the  procedure  of  offset  loading 
is  likewise  justified. 

The  final  results  of  the  tests  are  given  in  Figure  32  -  a  curve  of  critical 
load  versus  offset  distance.  On  this  curve  a  "theoretical  line"  has  been 
superposed,  derived  on  the  assumption  that  the  stresses  in  the  shell  can 
be  computed  from 


e  + « =  i  u  +  m 

A  +  Z  A  v  +  r / 


(2) 


and  from  the  hypothesis  that  f  =  f  =  a  constant  for  a  given  shell  when 
buckled.  The  agreement  betweenthe  observed  and  the  "predicted"  is  remarkable. 
Thus,  it  is  concluded  from  the  second  family  of  tests  that  a  cylindrical 
shell  under  combined  flexure  and  compression  buckles  when  the  compressive 
stress  reaches  a  level  which  would  have  initiated  buckling  under  pure  axial 
compression. 

This  result  is  in  agreement  with  the  result  of  Seide  and  Weingarten  but  in 
disagreement  with  that  of  Fliigge.  However,  it  is  only  fair  to  point  out 
that  Fliigge  qualifies  his  result  by  defining  a  buckle  aspect  ratio,  but  this 
ratio  does  not  appear  to  occur  in  practice. 

The  stress  obtained  as  the  critical  value  for  the  cylinder  under  pure  axial 
load  is  only  slightly  lower  than  that  given  by  the  classic  formula 


a  = 


(3) 


The  variation  is  not  sufficient  to  discredit  this  formula  or  to  justify 
Fischer's  value  of  O.065  x  O.605  Et/r.  It  is  the  normal  inaccuracy  due  to 
operator  and  equipment  error. 
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STRESS,  THOUSANDS  OF  psi 


Figure  23.  Stress-Strain  Curve  for  Material  of  the  Series  1 
Test  Specimen. 
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TABLE 

VI  -  continued 

Load 

Load 

Number  of 

Buckle 

Eccentricity,  e  Postition 

(lb) 

Buckles 

Increment 

3/8r  6 

462 

20 

482 

33 

13 

502 

53 

20 

522 

69 

16 

542 

80 

11 

7 

452 

16 

472 

4o 

24 

492 

51 

11 

512 

73 

22 

532 

84 

11 

8 

475 

25 

490 

38 

13 

505 

58 

20 

520 

72 

14 

535 

82 

10 

550 

89 

7 

l/2r  1 

400 

3 

420 

26 

23 

440 

52 

26 

460 

70 

18 

480 

76 

6 

2 

400 

7 

420 

23 

16 

440 

50 

27 

460 

70 

20 

480 

78 

8 

3 

385 

6 

405 

19 

13 

425 

46 

27 

445 

69 

23 

465 

77 

8 

4 

370 

6 

390 

24 

18 

410 

33 

19 

430 

50 

17 

450 

6l 

11 

470 

69 

8 

5 

402 

14 

422 

40 

20 

442 

58 

18 

462 

74 

16 

472 

80 

6 

6 

400 

9 

420 

2? 

13 

44o 

48 

26 

460 

69 

21 

480 

76 

7 

TABLE 

VI  -  continued 

Load 

Eccentricity,  e 

Position 

Load 

(lb) 

Number  of 
Buckles 

Buckle 

Increment 

l/2r 

7 

395 

11 

415 

35 

24 

435 

52 

17 

455 

66 

14 

475 

74 

8 

8 

390 

4 

410 

11 

7 

430 

39 

28 

450 

64 

25 

470 

77 

13 

i - 1 

%  OF  BUCKLE  POPULATION  AT  LOAD 


GRAPHIC  REPRESENTATION  OF  THE  DATA  OF  TABLE  E  NORMAL 
PROBABILITY  COORDINATES. 


POSITION  I 

520  540  2 

480  500  520  540  3 

LOAD,  L  B 


460 

480 

500 

520 

540 

460 

480 

500 

440 

460 

Figure  24,  Results  of  Circular  Traverse  of  Loading  for  the  Shell 

Specimen  of  Series  2  Tests,  l.cccntricity  of  Load  =  3/8r. 
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%  OF  BUCKLE  POPULATION  AT  LOAD 


480 

800 

520 

540 

POSITION 

6 

440 

460 

480 

500 

520 

540 

580 

4 

440 

480 

480 

500 

520  540 

5 

LOAD,  LB 

Figure  24.  Continued. 
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440 

460 

480 

800 

820 

840 

POSITION 

460 

480 

800 

920  840 

LOAD,  LB 

Figure  24.  Continued. 


_J _ I 

r 

960  6 
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%  OF  BUCKLE  POPULATION  AT  LOAD 


7.  OF  BUCKLE  POPULATION  AT  LOAD 


t 

5 


LOAD,  LB 

Figure  25.  Continued. 


4 

480  5 
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OF  BUCKLE  POPULATION  AT  L 


380  400  420  440  460  480  POSITION  6 

400  420  440  460  480  7 

380  400  4  20  440  460  8 

LOAD,  LB 

Figure  25,  Continued. 


MODE  VALUE  OF  THE  BUCKLE  DISTRIBUTION,  LB 


>vm»rw**»IW 


BASIC  DATA  CONTAINED  IN  TABLE  VI 

DETERMINATION  OF  MODE  VALUES  FOR  6>  3/8r 
IS  ILLUSTRATED  IN  FIGURE  27 


\  A* 

'  mS 


ECCENTRICITY  3/8r  MEAN  VALUE  «  490  LB 


’ - - - - - 0 - 

ECCENTRICITY 

l/2r  MEAN 

- - - 0 - ' 

VALUE >425  LB 

_ Q _ _ _ n _ £ 

- W - v - 

MODE  VALUE 

POSITION 

e=3/8r 

e=|/2r 

1 

490 

424  LB 

2 

492 

423 

3 

485 

420 

4 

492 

420 

5 

491 

424 

6 

494 

433 

7 

484 

426 

8 

_ 1 _ 1 _ 1 _ 

494 

_J _ L 

431 

_ 1 _ 1 _ 

LOAD  POSITION  DURING  CIRCULAR  TRAVERSE 


Figure  26.  Maximum  Buckling  Rate  Load  for  Liccentric  Loading  at  Light 
Equally  Spaced  Circular  Positions. 
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STANDARD  DEVIATION  OF  THE  BUCKLE  DISTRIBUTION, LB 


BASIC  DATA  CONTAINED  IN  TABLE  H 


0  I  234  S  6  7  8 


LOAD  POSITION  DURING  CIRCULAR  TRAVERSE 

Figure  27.  Initial  Buckling  Region  of  Shell  Specimen  of  Series  2 
Tests  as  Revealed  by  Circular  Traverses  of  Loading. 
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TABLE  VII.  CENTRAL  AXIS  LOADING  FOR  THE  SHELL  SPECIMEN  OF 
SERIES  2  TESTS 


Buckle  Count  vs  Load 


rOVO  CO  l/\CO  t— 
CIMOI  W  H  H 
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.... 


TABU:  VIII 

-  continued 

Load 

Eccentricity,  e 

Load 

(lb) 

Number  of 
Buckles 

Buckle 

Increments 

5/8r 

340 

4 

360 

13 

9 

380 

28 

15 

400 

43 

15 

420 

58 

15 

440 

69 

11 

3/4r 

305 

8 

325 

13 

5 

350 

34 

21 

370 

52 

18 

390 

64 

12 

410 

70 

6 
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%  OF  BUCKLE  POPULATION  AT  LOAD 


j 


500 

520 

540 

560 

3/8  r 

400 

420 

440 

460 

480 

l/2r 

360 

380 

400 

420 

440 

5/8r 

300 

320 

340 

360 

380 

ECCENTRIC  LOAD,  LB 
I'igure  28.  Continued. 


6l 


*aj 


r*  • 

« 

ry  »  . 
jd^*  ♦"'•  r* 

TOTAL  NUMBER  OF  BUCKLES 


LOAD  ECCENTRICITY 


Figure  31.  Total  Buckle  Population  as  a  Function  of  Load  Eccentricity, 
Shell  Specimen  of  Series  2  Tests, 
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MAXIMUM  BUCKLE  RATE 
LOAD,  LB 


0 


LOAD  ECCENTRICITY 


Figure  32.  Maximum  Buckle  Hate  Load  Versus  Load  Fcccntricity,  Shell 
Specimen  of  Series  2  Test*. 
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CONCLUSIONS 


The  work  reported  herein  is  substantial  experimented  evidence  that  unstiffened 
right  circular  thin-walled  cylindrical  shells  buckle  under  nonuniform  axled 
load  conditions  when  the  maximum  compressive  stress  in  the  shell  reaches 
the  level  which  would  cause  instability  under  uniform  load  conditions. 

It  is  demonstrated  that  sound  statistical  data  with  regard  to  shell  stability 
can  be  as  reliably  obtained  from  well-planned  tests  on  single  specimens  as 
from  a  multitude  of  teats  on  a  wide  range  of  specimens. 

A  strong  indication  is  given  that  there  is  no  need  to  use  elaborate  methods 
of  fabrication,  since  incomplete  shells  may  well  be  as  revealing  as  complete 
ones. 
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AFI'JNDxX 

STATISTICAL  TREATMENT  OF  THE  RESULTS  OF  THE  FIRST  SERIES  OF  TESTS 


The  data  obtained  from  buckling  tests  of  3^9  cylindrical  shell  specimens  are 
contained  in  Tables  I  through  V.  These  tests  are  described  in  detail  in  the 
body  of  the  report,  and  the  test  data  are  interpreted  graphically  in  Figures 
19  and  20.  The  data  of  Figure  19  pertain  to  loading  distributions  designated 
Type  A,  while  those  of  Figure  20  represent  the  results  from  loading  distri¬ 
butions  designated  Type  B.  The  purpose  of  this  appendix  is  to  present  the 
results  of  a  statistical  analysis  which  was  performed  on  the  data  of  these 
two  figures. 

An  underlying  physical  relationship  relating  the  fraction  of  cylinder  end 
perimeter  loaded  to  a  corresponding  value  of  buckling  load  is  unknown.  It  is 
true  that  the  work  of  Bijlaard  and  Gallagher,  discussed  in  the  Introduction  to 
this  report  and  described  in  detail  in  Reference  5,  provides  a  theoretical 
basis  for  predicting  such  a  physical  relationship.  However,  valid  theoretical 
work  covering  the  range  of  loading  distributions  of  interest  in  this  study 
could  not  be  found.  Thus,  the  statistical  problem  before  us  is  to  determine 
in  terms  of  confidence  levels  the  degree  of  association  which  exists  between 
the  fraction  of  perimeter  loaded  and  the  buckling  load. 

The  procedures  followed  are  those  outlined  in  the  text  by  Bowker  and  Lieber- 
man.  3  Similar  procedures  are  available  in  most  modern  references  on  experi¬ 
mental  statistics,  a  notable  example  being  Handbook  91  of  the  National  Bureau 
of  Standards.  A  confidence  level  of  95  percent  was  chosen  in  all  the  work 
presented  here. 

BUCKLING  DATA  FROM  TYPE  A  LOADING 


These  are  the  results  shewn  graphically  in  Figure  19.  The  general  appear¬ 
ance  of  the  plot  irxiicates  a  high  degree  of  linear  association  between  the 
fraction  of  perimeter  loaded  and  the  buckling  load.  Also,  it  appears  that 
a  single  line  representation  of  the  data  would  pass  through  the  origin. 

It  is  not  surprising  that  the  latter  tendency  should  appear.  An  examination 
of  details  of  the  Type  A  load  distribution  described  in  Figures  10,  11,  and 
12  shows  that  the  magnitude  of  load  required  to  produce  a  given  level  of 
stress  in  the  shell  wall  approaches  zero  as  the  fraction  of  perimeter  loaded 
approaches  zero. 

All  numerical  quantities  from  the  experimental  data  required  in  the  analyses 
presented  here  are  found  in  Table  IX. 

Test  for  Linearity 

At  four  values  for  the  fraction  of  perimeter  loaded,  more  than  one  mean 
value  of  buckling  load  was  observed.  Thus,  we  may  analyze  the  data  for 
linearity  using  an  F  test. 

Hypothesis:  the  association  between  the  variables  is  linear. 

Rejection  Criterion: 
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F  = 


n 

r 


k  l  <pi  -  pi>‘ 
1=1 _ 

n  -  2 


n 

T 

L 


v 

> 


(pij  -  pi> 


i=i  j=i 

“5fk  -  iT 


(4) 


^  ^  /,  , \  =  the  100a  percentage  point 

a;n  -  2,  n(k  -  1)  of  ^  p  ^lstrlb^loS  vlth 

degrees  of  freedom  n  -  2 
and  n(k  -  l).  0£  =  .05  at 
95 "percent  confidence  level. 


When  values  are  taken  from  Table  IX  to  compute  F, 

f  =  2.  1^63.  +  bj^oo  _  >922 


(5) 


F.05;  6,8  ~ 


3.58 


.922  <  3.58  and  the  hypothesis  of  linearity  is  not  rejected  at 
the  95 -percent  level. 

Test  for  Zero  Intercept 


Having  accepted  the  hypothesis  of  linear  association,  we  now  proceed  to  fit 
a  straight  line  to  the  data  points.  For  this  purpose  we  use  the  method  of 
least  squares. 

P  =  a  +  bf  a  =  vertical  <*xis  intercept  (6) 

b  =  slope 

n  n 

b  =  y  (fi  -  v&f  p)  +  y.  (fi  ■ 1)2 

i=l  i=l 


a  =  P  -  bf  =  1058  -  1780  x  .5^5  =  88  lbs 

Thus,  P  =  88  +  1780f  is  a  least-squares  linear  representation  of  the  assoc¬ 
iation  between  the  variables. 


In  view  of  the  physica]  requirement  described  above,  that  such  a  linear 
representation  must  extend  through  the  origin  as  a  limit,  it  would  appear  that 
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the  value  of  88  pounds  indicated  experimentally  for  the  vertical  axis 
intercept  is  accidental.  We  now  test  for  this  occurrence  using  a  Student 
"t"  test. 

Hypothesis:  a  =  0 
Rejection  Criterion: 


a 


*  *01/2;  n-2 


(7) 


where 


ta/2-n-°  ~  percentage  point  of  the  "t" 

'  ’  ”*■  distribution  with  degrees  of  freedom 
n-2 


and  the  hypothesis  that  a  =  0  is  not  rejected  at  the  95-percent  level. 
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Straight- Line  Fit  With  Zero  Intercept 


Accepting  the  hypothesis  of  linearity  and  zero  intercept  on  the  vertical  axis, 
we  now  fit  a  least-squares  line  through  the  experimental  points  which  extend 
through  the  origin. 


P  =  bf 


(8) 


b  = 


n 


i=l 


n 


i=l 


J899 

3.1010 


=  1890 


F  =  l890f 


This  is  the  line  shown  in  Figure  19  and  is  considered  to  be  the  best  statistical 
representation  of  these  experimental  data. 


Confidence  Interval  Estimate  of  Slope 


As  a  final  analysis  of  the  data  of  Figure  10,  a  95 -percent  confidence  interval 
estimate  for  the  slope  b  computed  in  Equation  (  8  )  will  be  made.  We  assume 
no  a/  priori  knowledge  of  values  for  standard  deviation  of  the  values  of  P  ; 
therefore,  we  base  the  estimate  on  a  Student  "t"  distribution.  ^ 


Interval  =  b 


±  t. 


'<y2;n-] 


(9) 


^f 


1  <fi  • 

i=l _ 

n  -  1 


49  lb 


/  )  f ?  =  ✓3.1010  =  1.76 

V  1.1 

\t  2;n  -  1  =  t.025;7  =  2'365 
Interval  =  1890  ±2.365  x 
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f 


=  1890  ±  66  lb 
[  102 1956  lb] 


Summarizing  Statement 

As  a  result  of  the  foregoing  treatment  of  the  data  of  Figure  19,  the  follow¬ 
ing  statement  can  be  made.  With  probability  .95,  the  true  experimental 
relationship  between  fraction  of  perimeter  loaded  and  buckling  load  for  Type 
A  loading  is  contained  within  the  bounds  of  the  two  lines  P  =  l82kf  and 
P  =  195 6f. 

BUCKLING  DATA  FROM  TYPE  B  LOADING 

These  results  are  shown  graphically  in  Figure  20.  As  was  the  case  for  Type  A 
loading,  a  high  degree  of  linear  association  between  fraction  of  perimeter 
loaded  and  buckling  load  is  indicated  by  the  plot.  However,  in  this  case  it 
is  clear  that  a  single  line  representation  of  the  data  would  not  pass  through 
the  origin.  Rather,  for  zero  fraction  of  perimeter  loaded,  a  relatively 
high  value  of  loading  is  apparently  required  to  cause  buckling. 

To  explain  the  latter  observation,  we  again  appeal  to  an  examination  of  the 
physical  details  of  how  an  external  load  is  transmitted  to  the  edge  of  the 
shell.  Figures  13  and  lk  show  the  Type  B  loading,  while  the  manner  in  which 
the  continuous  rims  participate  in  load  transfer  to  the  shell  is  described 
in  the  body  of  the  report.  It  is  clear  for  Type  B  loading  that  as  the  frac¬ 
tion  of  loaded  rims  perimeter  approaches  zero,  the  fraction  of  shell  perimeter 
reacting  this  load  approaches  a  limiting  value  greater  than  zero.  This  is 
due  tc  the  bending  stiffness  of  the  rims. 

In  the  case  of  a  concentrated  load  on  the  end  rim,  representing  zero  fraction 
of  perimeter  loaded,  the  shell  load  is  distributed  along  a  finite  length  of 
the  perimeter.  An  estimate  of  this  limiting  case  was  computed  by  treating 
the  end  rim  as  a  beam  on  an  elastic  foundation.  This  computed  value  was  then 
used  to  estimate  the  vertical  axis  intercepts  shown  in  Figure  20. 

All  numerical  quantities  from  the  experimental  data  required  in  the  follow¬ 
ing  analyses  are  found  in  Table  X. 

Straight-Line  Fit 

Acceptance  of  linearity  for  the  association  between  variables  is  based  on 
the  analysis  performed  for  Type  A  loading,  together  with  the  general  appear¬ 
ance  of  the  plot  in  Figure  20. 

The  rim  analysis  described  above  indicates  that  for  values  of  f  greater 
than  .15,  the  degree  to  which  rim  stiffness  participates  in  load  transfer 
to  the  shell  rim  diminishes  rapidly.  Thus,  a  straight-line  fit  will  be 
regarded  as  significant  only  for  values  of  f  greater  than  .15 • 
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I  TABLE  X.  TREATMENT  OF  EXPERIMENTAL  DATA  GIVEN 

IN  FIGURE  20 

"1 

i  i 

i 

fi 

Pi 

(P^P) 

1  (p,-*)  1 

MtX'iV 

1 

.15 

590 

-.40 

-620 

.1600 

248 

570 

20 

400 

2 

.25 

740 

-.30 

-470 

.0900 

l4l 

730 

10 

100 

3 

.33 

930 

-.22 

-280 

.0484 

62 

863 

67 

4500 

4 

.49 

1070 

-.06 

-l4o 

.0036 

8 

1115 

-45 

2030 

5 

.50 

1110 

-.05 

-100 

.0025 

5 

1130 

-20 

400 

6 

.67 

1300 

-.12 

90 

.0144 

11 

1396 

-96 

9220 

7 

•71 

1590 

-.16 

380 

.0256 

61 

1467 

123 

15120 

8 

.85 

1690 

-.30 

480 

.0900 

144 

1690 

0 

0 

9 

1.00 

i860 

-.45 

750 

.2025 

338 

1930 

70 

4900 

r 

L 

4.95 

10880 

.6370 

1018 

36670 

n  =  9 

9 

f  =  -  V  f  = 

n  L  1 

4.95 

9 

=  .550 

i=l 


9 

=  -  y  P  =  =  1210  lb  P  =  330  +  l600f 

,^1  y  (see  Straight-Line  Fit  section) 


f  =  Independent  variable  designating  fraction  of  perimeter  loaded, 
dimensionless 

P  =  Dependent  variable  designating  buckling  load,  pounds  of  fo^ge 
Pj=  The  value  of  buckling  load  observed  to  correspond  to  the  i 
value  of  f 


7k 


K 


n  n 

b  =  I  (fi  '  f)(Pi  -  *  L  (fi  ' 

i=l  i=l 

=  7§f  =  1600  lb 

a  =  P  -  bf  =  1210  -  1600  x  .550 
=  330  lb 
P  =  330  +  l600f 
Confidence  Intervals  for  P 


A  confidence  interval  estimate  for  values  of  P  for  a  given  value 
be  computed.  Again,  the  estimate  will  be  based  on  the  Student  " 
but ion. 


Interval  =  P  ±  t 


^72^-2  =  t.025;7 


2.365 


=  i  1  +  1.1+13(F  -  .55)2 


(10) 


of  f  will 
distri- 


(11) 
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f  >  .15 


(12) 


Interval  «  P  ±  57 >/l  +  1.4l3(f  -  .55)2, 

at  the  95 -percent  confidence  level. 

Three  sample  values: 

f  P  57  n/i  +  l.M3(f-.^)2  ,  Interval 

35  5  70  53  1b  [.507,633] 

.55  1210  57  [1153,1267] 

1.00  1930  65  [1865,19953 

Summarizing  Statement 

As  a  result  of  the  above  treatment  of  the  data  of  Figure  20,  the  following 
statement  can  be  made  for  values  of  f  greater  than  .15.  With  probability 
.95,  the  true  experimental  relationship  between  f  raction  of  perimeter  loaded 
and  buckling  load  for  Type  B  loading  is  contained  within  the  interval 

P  ±  57  /l  +  1.4l3(f  -  -55)2, 

P  «=  330  +  l600f 


f  >  .15 


(13) 


